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Preface 


There  are  two  objectives  I  want  to  bring  out  In  this  thesis. 

First,  the  main  objective  Is  to  verify  the  feasibility  of  the  new 
flight  testing  technique  for  the  Space  Shuttle  Orbiter  in  a  controlled 
environment.  This  Is  discussed  fully  In  the  last  section  of  this 
report.  Secondly,  I  attempt  to  summarize  the  data  reduction  method 
used  In  the  rest  of  the  report.  The  mathematical  formulation  may  be 
hard  to  follow,  especially  the  stochastic  estimation  section  of  this 
report.  I  recommend  that  the  reader  refer  to  Maybeck's  books  on  the 
subject.  I  found  It  easier  to  read  than  most  other  reference  material 
that  I  have  come  across.  Hopefully,  this  will  help  as  far  as  this 
section  Is  concerned,  and  the  rest  Is  straightforward. 

I  wish  to  express  my  gratitude  to  the  Air  Force  Wright  Aeronautical 
Laboratories  (AFWAL)  for  making  available  the  wind  tunnel  data  and  also 
the  descriptions  of  the  wind  tunnel  test  articles.  Without  the  material 
provided,  I  would  not  have  been  able  to  put  together  this  thesis.  I 
am  Indebted  to  my  thesis  advisor,  Capt  James  K.  Hodge  of  the  Department 
of  Aeronautics  and  Astronautics  at  the  Air  Force  Institute  of  Technology 
(AFIT).  His  helpful  suggestions  and  his  constant  encouragement  and 
understanding  has  made  all  this  possible.  I  also  would  like  to  thank 
Dr.  Will  Hankey  of  the  AFWAL  and  Capt  David  R.  Audley  of  the  Mathematics 
Department  at  AFIT  for  the  input  and  suggestions  that  kept  me  on  the 
right  track  throughout  this  thesis.  Last,  I  would  like  to  thank  Mrs. 

Cindy  McDonald  for  her  effort  in  typing  my  thesis. 
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Abstract 

A  transient  maneuver  testing  technique  has  been  developed  by 
the  Air  Force  Flight  Test  Center  (AFFTC)  for  aerothermodynamic  evalua¬ 
tion  of  the  Space  Shuttle  Orbiter  Thermal  Protection  System  (TPS) 
during  reentry.  The  objectives  are  to  determine  the  feasibility  of  the 
testing  technique  and  the  feasibility  of  using  the  transient  maneuver 
technique  for  wind  tunnel  heat  transfer  measurements.  The  transient 
maneuver  technique  can  enhance  the  capability  of  wind  tunnel  testing  by 
reducing  the  overall  required  testing  time,  and  in  turn  lower  the  cost 
of  testing  significantly. 

Wind  tunnel  data  were  obtained  using  the  TPS  materials  of  the 
Orbiter  as  wind  tunnel  test  articles.  Transient  and  steady  state 
deflection  angle  data  were  obtained  and  analyzed.  The  thermocouple 
measurements  along  with  the  wind  tunnel  conditions  were  used  as  inputs 
to  a  heat  estimation  computer  program.  The  program  named  "NEATEST"  was 
developed  by  AFFTC  and  used  system  identification  theory.  The  results 
of  the  estimates  were  compared  with  the  available  thin  skin  data  and  also 
with  the  Eckert  flat  plate  theory. 

The  principal  finding  indicates  the  heating  estimates  are 
significantly  lower  for  the  TPS  test  articles  than  the  thin  skin  test 
article  and  Eckert  theory.  The  discrepancy  is  believed  to  be  caused 
mainly  by  the  non- isothermal  wall  effect.  Finally,  the  transient 
maneuver  testing  technique  has  been  proven  since  the  transient  maneuvers 
agree  well  with  the  steady  state  results. 


I.  Introduction 


Background 

Today's  space  transportation  system,  the  Space  Shuttle  Orbiter, 
is  very  versatile.  One  unique  feature  about  the  Shuttle  is  its  reusable 
Thermal  Protection  System  (TPS),  which  consists  primarily  of  low  density 
ceramic  tiles,  (fief.  13)  For  example,  the  TPS  found  on  the  lower  surface 
of  the  Orbiter  are  High-Temperature,  Reusable  Surface  Insulation  (HRSl) 
tiles  made  of  a  special  kind  of  low  density  ceramic.  Flexible  Reusable 
Surface  Insulation  (FRSI)  made  of  coated  Nomex  felt  can  be  found  on  the 
upper  surface  where  temperatures  are  not  as  l. Although  the  TPS  has 
been  evaluated  by  ground  testing  in  a  controlled  environment,  it  was 
not  tested  in  a  real  flight  environment  before  the  first  flight  in 
April,  1981.  For  safety  reasons,  the  mission  of  the  Space  Shuttle 
Orbiter  has  been  limited  to  a  narrow  range  of  operations.  Therefore, 
a  primary  goal  is  to  expand  the  operational  envelope  of  the  Space  Shuttle. 

National  Aeronautics  and  Space  Administration  (NASA)  initiated  a 
flight  test  program  to  assess  the  capability  of  the  Space  Shuttle  Orbiter. 
The  first  several  flights  were  Initially  designed  for  flight  test  pur¬ 
poses.  The  first  flight  test  demonstrated  the  safe  operation  of  the 
Space  Shuttle.  The  subsequent  flights  are  carefully  designed  to  gradually 
expand  the  operational  envelope  of  the  Space  Shuttle  Orbiter.  Although 
the  assessment  program  consists  of  ascent,  orbital  and  reentry  phaset 
of  operations,  the  main  emphasis  in  this  report  will  be  on  the  reenti > 
phase.  (Ref.  6) 


with  the  baseline  assessment  program  in  mind,  a  new  flight  test 
technique  has  been  developed  by  the  Air  Force  Flight  Test  Center  (AfFTC) 
to  aid  In  the  evaluation  of  the  reentry  aerothermodynamic  environment  of 
the  Orbiter.  The  new  assessment  technique  consists  of  three  basic  steps. 
First,  wind  tunnel  data  and  a  flight  simulator  were  utilized  for  pre¬ 
liminary  flight  test  evaluation  and  planning.  Secondly,  a  new  data  acqui 
tion  technique  was  developed  which  combines  the  use  of  thermocouples  and 
a  transient  maneuver  testing  technique  to  obtain  flight  data.  Lastly,  a 
flight  data  reduction  technique  was  developed  which  utilizes  thermocouple 
measurements  to  estimate  the  aerothermodynamic  heating  of  the  Orbiter. 

The  overall  results  will  ensure  a  rapid  but  safe  envelope  expansion  of 
the  Space  Shuttle  Orbiter.  (Ref.  2,10,11) 

FI Ight  Test  Technique 

Mathematical  models  were  developed  by  AFFTC  for  the  flight 
simulator  at  AFFTC.  (Ref.  10)  Mathematical  models  were  developed  for 
seven  locations  on  the  Orbiter.  These  locations  are  ctitical  for  over¬ 
all  mission  evaluation.  They  are  assumed  Independent  from  one  another. 

In  addition,  results  at  one  point  are  not  extrapolated  to  any  other 
point  on  the  Orbiter  by  these  models.  The  locations  are  shown  in  Figure 
1.  Thermocouples  were  Imbedded  In  the  TPS  at  numerous  locations  on  the 
Orbiter.  (Ref.  17)  The  models  will  be  used  on  the  flight  simulator  to 
simulate  the  aerothermodynamic  environment  at  each  one  of  these  locations 
or  at  any  other  locations  which  may  be  critical.  The  results  from  the 
simulator  can  then  be  used  for  evaluation  and  planning  of  the  real  flight 
test  maneuvers.  In  addition,  the  mathematical  models  will  be  used  tor 
flight  data  analysis  purposes.  Similar  models  were  used  in  wind  tunnel 
data  analysis. 


Figure  1  Space  Shuttle  Control  Point  Locations 


Transient  flight  test  maneuvers  havv  been  and  will  be  performed 
during  the  Orbiter  test  flights,  and  flight  test  data  will  be  obtained 
under  these  conditions.  One  example  Is  the  pushover-pul  1  up  maneuver 
performed  during  the  Orbiter's  second  test  flight.  (Ref.  10)  One 
requirement  Is  that  the  duration  of  the  transient  maneuvers  should  be 
short  enough  to  not  affect  the  overall  trajectory  of  the  mission.  More¬ 
over,  the  maneuvers  are  being  gradually  expanded  on  subsequent  flights. 

The  thermocouple  data  from  the  flight  test  are  being  analyzed, 
and  a  digital  computer  Is  suited  for  this  purpose.  A  heat-estimate 
computer  program  called  “HEATEST"  has  been  developed  by  the  AFFTC.  The 
“HEATEST"  computer  program  utilizes  the  state  estimation  and  system 
Identification  theory  to  estimate  the  aerothermodynamic  parameters. 

The  system  Identification  theory  in  this  case  simply  uses  tlie 
available  Information  given  about  the  Input  and  the  output  of  the  pro¬ 
cess  and  uses  this  Information  to  estimate  the  parameters  of  interest 
In  the  process.  For  example,  in  Figure  2,  the  u  and  z  is  the  input  and 

V 

the  output  of  the  system  which  can  be  used  to  characterize  the  complete 
process  In  the  box. 


Figure  2.  A  Simple  Control  Input/Output  System  Model 

The  data  analysis  technique  consisted  of  a  state  predictor-cor¬ 
rector  and  a  parameter  estimation  scheme.  The  predictor-corrector  scheme 
requires  a  stochastic  process  model  which  was  developed  for  the  siinulator 


to  propagate  or  predict  the  temperatures  and  use  the  Kalman  estimatot 
or  filter  to  update  or  correct  the  temperature  of  the  model.  (Ref.  12,  14) 


The  parameters  are  estimated  by  maximizing  a  likelihood 
function.  (Ref.  1,  2,  8,  14,  16)  The  results  will  be  a  best  estin.ate 
of  the  parameters. 

At  last,  the  flight  testing  technique  combines  the  wind  tunnel 
and  the  flight- test  results  to  generate  an  enriched  simulator  data  base 
for  understanding  the  heating  behavior  of  the  Orbiter  during  reenii-y. 

Purpose 

This  flight  testing  technique  can  a  ic  ba  used  in  a  wind  tunnel. 

The  transient  maneuver  testing  technique  can  reduce  the  testing  time 
in  the  wind  tunnel  significantly.  This  will,  in  a  way,  enhance  the 
capability  of  wind  tunnel  testing.  The  purpose  of  this  thesis  is  to 
use  the  data  reduction  technique  which  was  developed  for  the  flight 
data  analysis,  and  the  wind  tunnel  experimental  data  obtained  by  the 
Air  Force  Wright  Aeronautical  Laboratories  (AFWAL)  to  verify  the  feasi¬ 
bility  of  the  transient  maneuvers  in  the  wind  tunnel.  Initially,  the 
wind  tunnel  test  was  conducted  to  verify  the  flight  testing  technique 
before  the  first  flight  test  maneuvers  were  committed.  The  reason  was 
that  the  wind  tunnel  data  could  be  obtained  in  a  controlled  envi roninent . 

For  example,  wind  tunnel  data  will  not  have  real  gas  effects,  limited 
data  sampling  rate,  eight-bit  noise,  and  the  time  skews  experitncc-ii 
in  the  flight  measurements.  For  this  reason,  the  results  of  tlie  wind 
tunnel  tests  could  answer  some  of  the  questions  about  the  fliqrii  ic (.tmiques 
that  were  not  available  before.  Based  on  these,  the  tlight  tesciii  i  tech¬ 
nique  can  then  be  improved  in  the  future. 


Overview 


In  the  following  sections,  the  mathematical  models  and  the 
stochastic  modelling  process  will  be  presented.  A  state  estimator  and 
a  parameter  estimation  technique  will  be  discussed  in  full  detail.  Th 
the  wind  tunnel  results  will  be  discussed. 


II.  Hathematical  Modelling  Process 


Introduction 

Hathematical  models  were  developed  for  a  flight  simulator  and  also 
used  for  flight  data  analysis.  In  this  case,  the  models  were  used  for 
wind  tunnel  data  analysis.  Models  for  the  aerodynamic  forced  convective 
heating  at  the  surface  of  the  TPS  and  the  thermal  conduction  through 
the  tiles  will  be  discussed.  These  two  models  will  be  treated  separately 
at  the  beginning  of  the  modelling  process.  In  the  end,  a  complete  state 
space  model  of  the  process  will  be  discussed. 

Heating  Model 

The  aerodynamic  heating  rate  at  the  surface  was  modelled  for  two 
TPS  test  articles  in  the  wind  tunnel.  The  convective  heat-transfer  can 
be  written  as 

q«h{TrUo-Ui)  2.1& 

or  in  terms  of  the  heat-transfer  film  coefficient, 

h  =  q/(RUo-Ui)  ^Mb 

where  is  the  stagnation  temperature,  Ui  is  the  surface  temperature, 
and  R  is  the  recovery  factor  which  is  assumed  to  be  0.9.  The  heating 
models  were  assumed  to  be  functions  of  the  f'^eestream  conditions  in  Lhc- 
wind  tunnel  test  section.  In  addition,  the  aerodynamic  heating  rate  was 
assumed  to  be  linear  in  deflection  angle.  This  is  not  a  bad  assumption 
based  on  previous  experiments  and  also  from  Eckert  flat  plate  theory. 
Therefore,  the  linearized  heating  model  can  be  formulated  as  follows; 


“f(a)  -  ho  +  ho(a-ao) 


where  h  Is  a  dimensionless  heat-transfer  film  coefficient  which  was 
non-dimensional ized  by  dividing  a  reference  heat-transfer  film  coefficient 
(Ref.  9,10,11,17) 


href  based  on.  Eckert  flat  plate  theory  for  laminar  flow  with  zero 
deflection  angle.  (Ref.  7)  Laminar  flow  is  appropriate  for  the  wind 
tunnel  test  which  was  conducted,  but  for  the  flight  data  analysis, 
turbulent  flow  might  be  used. 

In  equation  2.1,  h^  is  the  heating  derivative  with  respect  to 
the  angle  of  attack  a.  ho  is  the  zero  intercept  at  the  reference  value 

“0- 

The  mathematical  model  will  be  used  at  the  boundary  or  the  sur¬ 
face  condition  for  the  development  of  the  thermal  model.  Furthermore, 
the  heating  derivative  and  the  zero  intercept  are  the  primary  unknown 
to  be  estimated. 

Thermal  Model 

A  one-dimensional  conduction  model  was  assumed.  (Ref.  10,11) 

The  one-dimensional  assumption  was  suitable  for  two  reasons.  First, 
conduction  in  the  Space  Shuttle  TPS  is  very  low,  and  the  duration  of 
the  test  maneuvers  is  short  enough,  so  that  the  lateral  conduction  would 
be  small.  But  if  the  lateral  conduction  bec.ime  necessary  in  the  calcula¬ 
tion,  an  appropriate  scaling  factor  will  scale  the  thickness  of  ttic 
test  article  according  to  the  amount  of  conduction  generated  in  tlie 
lateral  direction.  (Ref.  10)  Second,  stochastic  estimation  was  usi a 
in  conjunction  with  the  measured  data  from  the  real  system  to  update 
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the  estimated  values.  Therefore,  the  error  generated  by  the  assumptiun 
will  be  eliminated. 

A  cross-section  of  the  thermal  model  is  shown  in  Figure  3.  For 
implementation  on  a  digital  computer,  n-discrete  node  points  were  est^b 
lished  to  approximate  the  temperature  profile  in  the  model.  There  are 
four  node  points  near  boundary  locations,  and  the  rest  are  distributed 
in  the  Reusable  Surface  Insulation  (RSI)  tile.  An  energy/heat  balance 
was  used  to  formulate  the  following  “n"  heat  transfer  equations 


CaPaAXaUi  =  -2(KA/aXA)Ui  +  2(Ka/aXa)L' 

-2oe(U^-U;i)  +  2f{a)href  2.3a 

(CaPaAXa  +  CbPbAX2)U2  =  2(Ka/aX^)Ui 

-2(Ka/aXa  +  Kb/aX2)U2  +  2(Kb/aX2)U3  2.3b 

CbPB(aX^.1  +  AX^)Ui  »  2(KB/AXi.i)Ui.j 

-2(KB/AXi.i  +  KB/AXi)Ui  +  2(KB/AXi+i 2. 3c 

(CbPbAXl_3  +  C(;PcAXj.)Ul.2  *  2(Kb/aXl.3)Ul_3 

-2(Kb/AXl.3  +  Kc/aXc)Ul.2  +  2(Kc/aXc)Ul.i  2.3d 

(CcPcAXc  +  CqPoAX[))Ul.j  “  2(Kq/AXc)Ul_2 

-2(Kc/AXc  +  +  2'iCD/AX[j)UL  2.je 

CdPqAXdUl  =  2{Kd/aXd)Ul.i  -2(K[)/AXp)UL  2.3f 


where  Cj,  Pj  and  Kj,  j  ■  A,B,C,D,  are  the  material  properties.  Ttie  material 
specific  heat  Cj  and  the  thermal  conductivity  Kj  are  a  function  of  tt.e  temp¬ 
erature  at  each  node  point,  and  the  thermal  conductivity  is  also  u  function 
of  the  local  pressure  in  the  tile.  Therefore,  Kj  is  apinox imated  tiie 
average  temperature  between  nodes,  and  Cj  is  approximated  at  each  node. 

is  the  temperature  at  each  node,  and  represents  u  time  rate  oi 
change  of  temperature.  Last,  aX,-  is  the  distance  between  nodes. 


The  primary  unknowns  are  mostly  in  Equation  3.2a,  the  surface 
boundary  condition.  The  unknowns  are  the  thermal  emissivity  e,  the 


aerodynamic  heating  rate  f(a).  and  the  coating  thickness  aX/^.  For 
example,  the  emissivity  can  be  as  low  as  .75  and  as  high  as  .92.  I  lie 
thermal  conductivity  Kj  was  originally  assumed  to  be  “a  priori"  knuwn 
parameter,  but  since  the  thermal  conductivity  depended  on  the  pressure 
inside  the  TPS'tile,  and  the  pressure  sensor  in  the  TPS  tile  might  not 
be  able  to  measure  the  pressure  correctly,  the  thermal  conductivi  ty 
is  also  treated  as  an  unknown.  The  specific  heat  Cj  and  the  density 
pj  are  assumed  to  be  a  priori  known  based  on  the  design  specification. 
Compl ete  State  Space  Model s  with  Stochastic  Process 

In  order  to  determine  the  unknown  parameters,  system  identitica 
tion  theory  was  used.  One  requirement  is  to  introduce  a  stochastic 
modelling  process.  The  model's  uncertainty  is  associated  with  the 
approximation  about  the  real  process.  The  n-di fferential  equations 
which  were  developed  earlier  can  be  arranged  in  a  more  convenient  state 
space  format.  One  reason  for  establishing  a  linear  state  space  model 
is  that  a  solution  of  the  state  can  be  easily  found. 

The  radiation  term  in  the  boundary  equation  is  highly  nonlinear. 
In  order  to  set  up  a  linear  state  space  form,  the  fourth  order  teni 
approximated  by  (Ref.  3,  18) 

Ulf(tn)  s  4{Ui(tn-i))3  Ui(tn)  -  3(011^.1))"  2.4 

where  Ui{t^_j)  is  a  known  temperature  which  is  propaga*^ed  from  Hit 
previous  time.  With  this  approximation,  the  n-di fferential  equuiiufi 


can  be  written  in  the  form 


U(t)  -  AU(t)  +  B  +  W(t)  ^.5 

where  A  is  the  n  x  n  tridiagonal  matrix  of  the  material  properties, 
and  ^  is  the  n-dimensional  vector  of  the  deterministic  inputs.  (Set- 
Appendix  A)  U(t)  is  an  n-dimensional  state  vector  at  time  t.  W(t) 
is  an  n-dimensional  white  Gaussian  noise  vector  with  statistics 

E{W(t)l  =0  2.6a 

E{W{t)W(T)^l  =  Q(t)6(t-T)  2.6b 

where  E{*}  is  the  expectation,  and  Equation  2.6  represents  a  zero  mean 
with  covariance  Q(t)6(t-T).  Q(t)  is  an  m  x  matrix,  and  6(t-T)  is  the 
Dirac  delta  function.  The  property  of  the  Dirac  delta  function  is  simply 
zero  everywhere  except  at  T*t.  (Ref.  2,5,12,14) 

There  are  reasons  why  white  Gaussian  noise  was  used.  White  noise 
is  uncorrelated  in  time,  and  has  equal  power  in  all  frequencies.  Because 
equal  power  in  all  frequencies  means  infinite  power,  white  noise  therefore 
can  not  exist  in  nature.  There  are  two  reasons  why  white  noise  can  be 
used.  First,  a  real  system  noise  especially  with  a  wideband  noise  usually 
can  be  approximated  quite  well  with  white  noise.  Second,  the  mathematics 
is  greatly  simplified  with  white  noise. 

A  Gaussian  process  will  ensure  a  "best"  estimate  since  the  mean, 
the  maximum  likelihood,  and  the  median  all  coincide.  In  addition,  only 
the  first  two  moments  are  necessary  to  completely  describe  the  statis¬ 
tics  in  the  process.  The  first  two  moments  are  the  mean  and  the  covariance 
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of  t+je  process.  Therefore,  with  a  white  Gaussian  process,  the  coi.iplixity 
of  the  process  can  be  considerably  reduced.  (Ref.  5.12.14) 

Initial  conditions  are  needed  to  completely  characterize  t(ie 
state  equations,  but  in  many  cases,  initial  conditions  may  not  be 
known  precisely.  Therefore,  the  initial  temperature  U(to)  are  charac- 


terized  by  the  statistics 

E{U(to)l  =  ^ 

/.7a 

E{CU(to)  -  ^][U(to)-^]^}  =  Pq 

/.7b 

where  ^  is  the  mean,  and  Pq  is  a  n  x  n  covariance  matrix.  The  initial 
noise  is  also  assumed  to  be  a  Gaussian  process.  In  addition,  the  be-ot 
estimate  of  £(tQ)  is  assumed  to  be 

To  complete  the  modelling  process,  a  measurement  model  is 
developed, 

Y{ti)  *  HU(ti)  +  V(t^)  ^.8 

where  Y^(ti)  is  k-vector  measurements  at  a  discrete  time  ti.  H  is  an 
k  x  n  matrix  which  locates  measurements  relative  to  other  nodes,  and 
V(t^)  is  a  "stationary"  white  Gaussian  process.  The  statistics  of  tiie 
k-vector  is 

E{V(tT)}=0  f'.9a 

T  ^  ^i  ” 

E{V(ti)v'(tj))  =  ^  ..9b 

~  ■  0  ti  tj 

where  R  is  a  time  invariant  k  x  k  matrix  for  stationary  white  Caussi.iii 
process.  All  sensors  are  assumed  to  be  independent  with  each  oilier. 
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In  summary,  combining  the  measurement  data  and  the  statistical 
Information,  a  best  estimate  of  the  state  car,  po  obtained.  Since  the 
measurement  data  Is  obtained  from  the  real  system,  unknown  parameters 
can  be  Identified  by  solving  the  maximum  likelihood  criterion  with  an 
iteration  scheme.  The  results  will  be  a  best  possible  estimate. 


III.  state  and  Parameter  Estimating  Technique 


Kalman  Estimator 

The  temperature  profile  can  be  estimated  using  a  Kalman  filter 
algorithm.  This  will  be  evident  in  the  next  section  when  solving  for 
the  maximum  likelihood  criterion.  In  order  to  estimate  the  temperature 
profile,  the  values  of  the  unknown  parameters  have  to  be  known,  but  ror 
now,  the  parameters  are  assumed  to  be  “a  priori"  known  values  in  which 


U*(t^)  =  U(t/) 


6  =0* 


j.l 


where  ^*(t^)  is  assumed  to  be  the  "optimar*  .  imate ,  and  L)(tj^)  is  the 
best  estimate  from  the  Kalman  estimator  using  the  a  priori  known  values 
of  the  unknown  parameters. 

The  estimator  Involved  a  predictor-corrector  scheme.  The  temp¬ 
erature  profile  at  time  t,  is  predicted  by  the  heat  and  thermal  models 
and  the  a  priori  Information  at  time  t^.^.  It  is  then  updated  by  the 
measurement  data  at  time  t^  which  is  provided  by  the  real  system.  The 
temperature  profile  is  propagated  with  the  same  scheme  through  the  entire 
time  history.  Therefore,  the  Kalman  filter  equations  can  be  formulated 
as  follows: 


U(ti’)  =  ♦(ti  -  ti.i)U(ti./)  + 


T )^(t )dT 


P(ti‘)  “♦(t^  -  t^,j)P(t^_j^)  i^ti  -  ti.l) 


•t)Q(T)4(t  j 


T  )d  I 


..  i 
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where  Equation  3.2  and  3.3  are  the  equivalent  difference  equation 
of  Equation  2.13  and  Equation  2.14.  ^(tj")  and  P(t^")  are  the  predicted 
temperature  profile  and  the  covariance  at  time  t^  using  the  updutod  values 
U.{ti_i  )  and  )  from  the  previous  time  t^.j. 

The  updated  equations  Involve  solving  the  conditional  density 
function 

conditioned  on  the  entire  measurement  sequence  Y^(ti)  up  to  time  t^.  Ihe 
derivation  can  be  found  in  Maybeck,  and  the  final  results  can  be  Mriticn 
as  (Ref.  2) 

U(t/)  -  U(t^')  +  K(ti)[lJi  -  HU(t.")]  J.i; 

P(ti^)  =  [I  -  K(ti)H]P(ti")[I  -  K(ti)H]^ 

+  K(ti)P(ti‘)K(ti)^  3.6 

where  the  "Kalman  gain" 

K{t^)  =  P{t^")H^[HP(t^')H^  +  R]‘^  J.7 

U(t^  )  and  P(t^  )  are  the  "a  posteriori"  estimates  of  the  state  and  the 
covariance.  Since  the  derivation  of  Equation  3.2  to  Equation  3.6  have 
assumed  the  values  of  the  unknown  parameters  are  "apriori"  known  and 
which  is  not  the  case,  a  parameter  estimation  technique  using  system 
identification  theory  was  introduced. 

Maximum  Likelihood  Criterion 

There  are  many  methods  which  can  be  used  for  estimating  puraim-ters. 
A  maximum  likelihood  technique  was  chosen  in  this  case  for  several  n  asuns 
First,  in  certain  cases,  maximum  likelihood  technique  had  been  shown  to 


have  nice  asymptotic  properties.  For  example,  with  a  nice  probabiliiy 
density  function,  the  mean,  mode,  and  the  median  are  usually  very  close 
together,  and  with  these  conditions,  the  maximum  likelihood  estimate 
will  be  very  close  to  the  "best”  or  the  "optimal"  value.  Furthermot e , 
the  maximum  likelihood  method  can  be  applied  to  a  large  variety  of  cases, 
and  many  other  estimation  techniques  are  just  a  special  case  of  the 
maximum  likelihood  technique.  (Ref.  1,14,16) 

In  order  to  establish  the  maximum  likelihood  criterion,  the  first 
step  is  to  find  a  likelihood  function  that  will  be  best  suited  for  the 
estimation  purpose.  There  are  many  likelihood  functions  which  can  be 
used,  but  one  that  has  been  found  to  be  the  best  is  (Ref.  2,14). 

^U(t^),Y(ti)/0^i*^i^®^  3.8 

where  Equation  3.8  is  a  joint  probability  density  function  of  U(ti)  and 
l[(t^)  conditioned  on  0  .  (Ref.  2)  U(ti)  is  the  temperature  profile 
at  time  t'j.  Y(ti)  is  the  measurements  sequence  up  to  time  ti ,  and  9 
is  a  set  of  unknown  parameters. 

U(ti)  -  [Ui(ti)U2(t^)...Un(t.)]^ 

I(t^)  »  [Y(ti),Y(t2)  •••  Y(ti)]^ 

0  *  [AX^f  ^0»^a3^ 

{  is  a  set  of  acquired  measurement  values  to  be  used  as  data 

Hi  “  [yi .  y2  •••  yj] 

i  and  0  are  the  dummy  variables  corresponding  to  U(t-)  and  0. 


r-  - . 


m  * 


i  e 


r'.: 


One  way  of  maximizing  the  likelihood  function  Is  to  take  the 
gradients  of  Equation  3.8  with  respect  to  the  temperature  state  £, 
and  the  unknown  parameters  e.  This  is  the  same  as  finding  an  optimal 
set  of  ^(t^)  and  0  at  time  t'j  that  Is  most  likely  to  produce  the 
measurements  ^(ti).  The  gradients  of  the  likelihood  function  witli 
respect  to  K  and  e  can  be  written  as 

^T 


V^LU.e;  Hi] 


i  “ 

0  «  0* 


=  O' 


3.yd 


VQL[c.e;  Hi] 


£  “  U(ti)* 
6  -  0* 


=  O' 


3. '3b 


where  L[£,eiHi]  is  the  natural  logarithm  of  Equation  3.8, 

LU.e;  HO  =  ^'’^^U(t^),Y(t^)/0  (i»Hi/®)) 


3.10 


A  closed  form  solution  for  Equation  3.9  may  not  always  be  available. 
Therefore,  an  Iterative  algorithm  was  generated  to  solve  the  criterion 
presented  In  Equation  3.9.  The  iterative  algorithm  can  be  easily  imile- 
mented  on  a  digital  computer. 

Solution  of  the  Maximum  Likelihood  Criterion 

In  order  to  satisfy  the  maximum  likeln.ocd  criterion,  first 
Bayes*  rule  Is  applied  to  Equation  3.8  repeatedly  through  the  entire 
measurements  sequence  to  obtain 


I'J 


'U(t^),Y(t^)  =  'U(t^)/Y(t^)-'Y(t,.) 

=  ^U(t^)/Y{t.).^Y{t.)/Y; 


“  Vt^)/Y(t^)  n  ^Y(tj)/Y(tj_^)  3.11 

j-1 

where  6  has  been  suppressed  to  simplify  the  notations,  and  ^  at  the 
first  measurement  is  not  conditioned  on  any  previous  lueasureiiient .  Cach 
separate  density  in  Equation  3.11  can  be  written  as 

^U(t^)/Y(t.)  =  [2ii"''2p{t/)ir*exp(-i 

{•}  =  -Ui-  U(t/)]^P(t/)-'U  -  U(t/)]  J.12 

where  U(t^*)  and  P(t^^)  are  the  mean  and  covariance  which  are  conditioned 
on  the  measurement  sequence  _Y{t^.)  at  time  t^  and  the  unknown  parameters 
0,  and  similarly 

,)  ■  1  exp  {.1 

{•)=  -iCHj  -  HU(tj)]V(tj)-‘[lJ.  -  HU(t.)]  3.13 

where 

F(t.)  =  HP(t/)H^  +  R  y.14 

J  w 

By  substituting  Equation  3.12  and  Equation  3.13  back  into  Equation  j.tu, 
and  solving  for  Equation  3.9a,  the  solution  of  the  gradients  of  the  like¬ 
lihood  function  with  respect  to  ^  is 

-  U  -  0(t,*)fp(t  *)•'  -  0^ 

-  "  '  '  <  ■  u*(c)  " 


or 


u’U,)  •  u(t/) 


I  9  -e'd,)  3.15 

From  looking  at  Equation  3.15,  the  maximum  likelihood  estimate  of  tin- 
state  obtained  by  using  the  Kalman  filter  algorithm  pro¬ 

vided  that  the  values  of  the  parameters  0  are  the  maximum  likelihooJ 
estimates  of  the  parameters  0*. 

The  same  procedures  can  be  applied  to  Equation  3.9b,  the  gradient 
of  the  likelihood  functions  with  respect  to  6.  The  solution  of  Equation 
3.9b  is  in  Maybeck,  Vol .  II,  Chapter  10.  A  closed  form  solution  of 
Equation  3.9b  generally  is  not  possible.  An  iterative  algorithm  will 
have  to  be  used. 


Newton-Raphson  A1 gorithm 

The  Newton-Raphson  technique  was  used  to  find  the  best  estimates 
of  the  parameters  0*  denoted  by  0*{t^).  The  Newton-Raphson  equation 
can  be  written  as 


0*(t^)  -  0*(t^)  -  v2L[G*{t^).0*(t.);g.]"^ 

•  VeL[U*(t.).e*(t^);iYi^ 

^  A  as 

where  e*(t^)  is  the  previous  estimate,  and  L[U*(t^ ) ,0*(t^ ) ;  is  the 
likelihood  function  evaluated  using  U*(tp  (the  best  estimate  of  the 
temperature  profile  from  the  Kalman  estimator  using  the  best  available 

A  A 

parameter  estimates  eA(t^))and  u*(t^). 
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Unfortunately.  Equation  3.16  involves  solving  a  second  order 
partials  differential  equation,  and  the  computation  will  be  enormous. 
An  approximation  called  “scoring"  was  suggested  by  Rao  (Ref.  15)  which 


simplifies  the  computation  substantially.  The  accuracy  of  the  compulation 
can  be  maintained  as  long  as  the  "number  of  samples  is  large".  (Ref.  2,14,15) 
The  approximation  is 


V§L[U*(t.).6*(t^);a.]  =-J[t.,y*(t.),  o*(t.)  ;y.] 


3.17 


where  J[‘]  is  called  the  conditional  information  matrix  given  by 


J[t^.U*(t^).0*(t^);  *  E{VQLU*(t.),6;  g.] 

%L[U*(t.),9, 


O-u 


3.18 


Therefore,  with  the  approximation,  the  computation  of  Equation  3.16  will 
only  involve  the  evaluation  of  the  first  order  partials,  and  Equation 
(3.16)  becomes 

-1 


0*(t^.)  =  0*(t^)  +  O[t^.U*(t^.).0*(t^);l|p' 
•  VQL[U*{t^),0*(t^);y^3 


3.19 


where  the  term  of  the  first  order  derivative  can  be  approximaLud  by 
(Ref.  2) 

i 

_2_L[!!.(t().e.(t():  M)]  “  3U(tj  ■  {■) 

.  J-i  ■ 


^  T 

where  8y(tj“)  is  the  sensitivity  which  can  be  evaluated  by  taking 

aeT” 

the  partials  of  Equation  2.3  with  respect  to  the  parameter  (see 
Appendix  B)  where  K  *  1,2,3  •••  P. 

The  "conditional  information  matrix"  can  also  be  approxiinaied 
by 

Jk,  .  E  30(t,-)^  HV‘(tj)H  3U(tj  )  3  20 

30|^  301 

where  is  the  kl^^  element  of  the  P  x  P  matrix.  The  3y(tj") 

is  the  "a  priori"  sensitivity  propagated  from  "a  posteriori"  sensUivity 
ayCtj.!"*^)  or  simply 

aok 

3y(tj-)  »  3y(tj.i+)  ^ 

l0k  iek 

The  "a  posteriori"  sensitivity  can  be  evaluated  by 

3y(tj'^)  =  [I  -  K(tj)H]  3y(tj~)  3  22 

30k  00k 


Therefore,  the  sensitivity  is  updated  and  propagated  through  the  entire 
sample  time  history.  With  the  information  provided  by  the  "score"  vector 
and  the  "conditional  information  matrix",  the  best  estimate  of  the  para¬ 
meters  can  be  found.  Finally,  the  best  estimate  of  the  parameters  will  . 
then  be  used  in  the  Kalman  state  estimator  to  propagate  the  temiiei  aures. 
The  iteration  will  continue  until  the  last  1 1 .m  .i i  iuii  is  reachL-d,  jnd 
hopefully,  the  parameters  will  already  be  converged  to  the  optmui  values. 
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A  Summary  of  the  Data  Reduction  Technique 


The  whole  data  reduction  technique  Is  SL.marIzed  In  Figure  4. 

The  process  is  characterized  by  four  steps.  First,  initial  conditions 
usually  were  generated  by  first  running  a  small  segment  of  the  maneuver 
and  using  the  output  of  the  process  as  initial  conditions.  Constant 
initial  conditions  are  assumed  for  the  first  small  segment,  and 

A 

Initial  sensitivities  V-U(t„)  are  assumed  to  be  zero.  The  initial 
covariance  for  the  second  time  segment  are  usually  reinstated  to  the 
values  given  In  the  error  models.  Second,  with  these  initial  condi¬ 
tions,  temperatures,  covariance,  and  sensitivities  can  be  estimated 
using  Equation  3.2,  Equation  3.3,  and  the  sensitivities  equation  derived 
in  Appendix  B.  Since  Equation  3.2  only  involved  solving  the  deternh n i stic 
part  of  the  model,  a  method  described  in  Williams  (Ref.  18)  can  be  used. 
First,  the  time  rate  of  change  of  temperature  is  approximated  with  a  first 
order  backward  difference  given  by 

U  »  [U(tJ  -  U(t^_l)]/at  3.23 

Since  the  radiation  term  is  the  only  nonlinear  term,  a  special  iterative 
method  was  used.  A  Newton-Raphson  iteration  technique  with  a  quadratic 
extrapolation  scheme  which  Is  described  fully  in  Williams  was  used  to 
solve  Equation  2.3a,  the  surface  boundary  condition.  Therefore,  with 
the  surface  temperature  known,  the  rest  of  the  temperature  profile,  n-I, 
can  be  solved.  This  method  has  been  found  to  be  more  efficient  thaii 
solving  Equation  3.2  which  has  also  been  linearized. 


The  next  step  is  to  update  the  estimates  using  Equation  3.‘j  lu 
Equation  3.7  and  Equation  3.22.  This  procedure  is  repeated  and  propagated 
up  to  time  t^ ,  the  end  of  a  time  segment.  In  addition,  the  score  am!  the 
conditional  information  matrix  will  be  generated  by  adding  up  the  sensi¬ 
tivity  estimates.  Last,  the  parameters  can  then  be  updated  with  the 
information  provided  by  the  score  and  the  conditional  information  mattix. 
Finally,  the  whole  process  is  then  started  all  over  again  and  tliis  tine 
with  the  new  initial  conditions. 


IV.  Wind  Tunnel  Result 


Physical  Winfl  Tunnel  Test  At  tides  Descriptions 

A  wind  tunnel  test  was  conducted  by  the  AFWAL  in  support  ot 
AFFTC.  The  experiment  is  conducted  in  a  Mach  14  wind  tunnel,  r-!,- 
purpose  is  to  simulate  the  flow  field  of  the  Space  Shuttle  Orbitei  dur¬ 
ing  reentry. 

The  primary  wind  tunnel  model  consisted  of  a  15  inche  lorn, , 

7  inche  wide,  and  2  inche  thick  stainless  steel  plate.  One  inipi.rtdnt 
feature  in  this  model  is  a  6  inch  square  cut-out  which  is  used  to  ould 
test  samples.  The  test  articles  consisted  of  the  Space  Shuttle  Hi '.>1 
tile  and  the  FRSI  material  and  thin  skin  stainless  steel  article.  (See 
Figures  5,  6,  7)  Thermocouples  are  imbedded  at  numerous  locations  in 
the  test  articles  as  shown  in  the  figures.  The  complete  flat  plute 
model  is  then  mounted  on  a  support  bracket  which  in  turn  is  inount< d 
on  the  pitch  sector. 

A  remote  shock  generator  was  also  used  for  some  test  runs. 

Since  shock  generated  data  was  not  used,  a  description  of  this  model 
was  not  necessary.  A  water  cooled  surface  which  is  shown  in  rifjuin  h 
was  used  with  Space  Shuttle  articles.  This  article  is  mounted  above 
the  flat  plate  model.  This  surface  will  provide  a  known  temperature 
which  heat  from  the  flat  plate  model  is  radiated  to. 

Therefore,  with  the  primary  flat  plate  model,  thermocouple 
measurements  were  obtained  by  pitching  the  model  with  transient  ai.d 
steady  state  maneuvers. 


PRESSURE 


HRSI  Wind  Tunnel  Results  and  Discussion 


A  systematic  approach  was  used  to  reduce  the  HRSI  data.  Then 
are  three  steps  in  this  process.  First,  a  transient  maneuver  of  the 
HRSI  data  will  be  reduced  to  estimate  the  heating  and  thermal  parameters. 
Second,  the  steady  state  and  transient  data  will  be  reduced  and  compared. 
Last,  a  comparison  between  the  HRSI  data  and  the  available  thin  skin  data 
and  also  with  the  Eckert  theory  should  demonstrate  the  feasibility  ot  the 
transient  maneuver  technique. 

A  transient  maneuver  was  analyzed.  The  maneuver  time  history 
and  the  temperature  time  history  are  shown  in  Figure  D-1.  The  maneuv.  r 
pitched  the  model  from  14  degrees  to  3  degrees  continuously,  (See 
Appendix  C)  The  main  purpose  for  analyzing  this  maneuver  was  to  identify 
the  necessary  heating  parameters  and  the  thermal  parameters.  A  surface 
thermocouple  was  used  for  implementation  to  "HEATEST".  For  some  reason, 
the  parameters  were  trading  with  one  another,  and  the  "HEATEST"  wasn't 
able  to  estimate  the  parameters.  One  reason  might  be  that  error  is 
present  in  the  measurements  that  was  not  accounted  for  in  the  mathematical 
model s . 

One  way  to  show  if  error  exists  in  the  measurements  is  to  let 
the  temperature  propagate  without  updating  by  the  Kalman  filter.  Heat¬ 
ing  and  thermal  parameters  were  fixed  similar  to  a  simulator.  Several 
values  for  the  parameters  were  used  to  generate  the  temperature  time 
history.  The  temperature  time  history  that  tracked  closely  with  the 
thermocouple  measurements  is  shown  n>  tiguie  b-  ITu.  time  lustory 
shows  a  lag  in  the  thermocouple  measurements  of  appruximately  tui..e 
seconds.  Three  seconds  lag  is  unrealistically  high  to  be  accounted  t>.i 
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by  anything  in  the  present  models.  This  has  shown  that  errors  do  i  Aist 
in  the  system,  and  the  error  could  be  in  either  the  inathematicd I  inodels 
or  in  the  measurements  obtained  from  the  real  system.  Further  investi¬ 
gation  will  be  necessary. 

Small  time  segments  of  maneuver  one  were  (jenerated,  and  the  results 

were  correlated  manually.  The  segments  are  shown  in  Figure  D-3.  (Jin-  assui.tp 

tion  that  was  made  is  that  segments  containing  steady  state  deflection 
angle  should  correlate  well  with  one  another.  Furthermore,  segments 
that  are  away  from  the  wind  tunnel  start  up  and  away  from  the  transient 
maneuver  should  correlate  very  well.  This  assumption  came  from  inspec¬ 
tion  of  Figure  D-2.  The  time  history  tracks  very  well  with  tlie  measure¬ 
ment  data  except  at  the  initial  wind  tunnel  start  up  and  at  the  tcui sient 

maneuver.  The  results  of  the  analysis  are  shown  in  Figure  D-4  to  ! igure 

D-6.  Several  interesting  results  were  apparent.  First,  the  results 
indicate  that  the  coating  thickness  is  best  approximated  to  be  .OUi  tt 
since  further  decrease  of  the  thickness  will  not  change  the  correl  tion 
very  much.  The  thickness  was  expected  to  be  between  .000833  ft  and  .00125 
ft.  However,  this  test  article  has  been  tested  numerous  times  in  an  arc 
jet.  Although  the  segments  correlate  well  compared  with  the  other  con¬ 
figurations,  the  transient  error  still  shows  up  in  the  figures.  Si.cond, 
the  only  effect  that  shows  up  in  the  figures  by  changing  the  theniul  con¬ 
ductivity  factor  is  the  magnitude.  The  wo  .i  case  based  on  atmospneric 
pressure  as  opposed  to  the  pressure  measurement  at  the  base  of  the  tile  , 
seems  to  be  with  a  thermal  conductivity  factor  of  three.  This  also 
shows  up  in  the  steady  state  correlation  which  is  shown  in  Figure' 


and  Figure  D-8.  Therefore,  the  convective  heating  parameters  uerc  -Limdled 
with  coating  thickness  of  .001  ft  and  thermal  conductivity  factur  -t  one. 

The  analysis  so  far  indicates  that  the  error  is  generated  m  the 
measurement  data.  Conversely,  this  could  mean  that  the  real  systti.i  was 
not  modelled  correctly  by  the  present  mathematical  model. 

The  next  step  was  to  reduce  the  steady  state  runs.  The  time  history 
is  shown  in  Figure  D-9  and  Figure  D-10.  One  attempt  was  made  to  estimate 
the  coating  thickness.  The  estimates  were  in  the  range  of  .OOU97  it  and 
.00102  ft.  The  results  further  confirm  the  coating  thickness  estii.iate 
from  the  previous  investigation.  An  atten,  t  .uS  made  also  to  estimate 
the  thermal  conductivity  factor,  but  was  not  successful.  To  furtiier 
confirm  the  coating  thickness  estimate,  the  Kalman  filter  was  turiud 
to  see  if  the  optimal  estimate  was  indeed  reached.  The  result  of  ilie 
steady  state  run  with  14  degree  deflection  angle  is  tabulated  in  Table 
D- I.  The  results  indicate  that  the  coating  thickness  estimate  was  an 
optimal  estimate. 

Several  transient  maneuvers  were  also  analyzed  using  the  esiimated 
thermal  parameters.  The  time  history  of  the  maneuvers  are  sliown  n. 

Figure  0-11  to  Figure  0-13.  The  procedure  used  in  tiie  first  analysis 
with  the  heating  and  thermal  parameters  fixed  similar  to  a  simulatir 
run  was  also  used  to  analyze  the  maneuvers.  The  lag  in  the  tempeijture 
time  history  also  shows  up  in  all  the  maneuvers  analysis.  (See  l  i  ,vire 
0-14  to  Figure  0-15)  By  looking  at  these  figures,  the  lag  seemed  ' 
be  reduced  with  less  maneuvers.  Even  though  the  lag  didn't  show  a.,  in 
Figure  D-l5,  the  lag  is  there,  but  it  is  small.  This  should  cinii:. 
that  the  lag  in  the  measurements  is  caused  by  tin.  liansient  ett.ct 

J4 


The  results  of  the  transient  and  steady  state  were  compared  with 
the  available  thin  skin  data  and  also  with  the  Eckert  flat  plate  theory. 
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The  comparison  is  shown  in  Figure  D-16,  First  of  all,  the  comparison 
between  transient  and  steady  state  were  relatively  good.  This  indicates  that 
steady  state  testing  technique  can  be  replaced  by  the  transient  testing 
method.  But  the  comparison  between  HRSl  results  and  thin  skin  data 
and  Eckert  flat  plate  theory  indicates  a  large  discrepancy.  There  are 
two  possibilities  that  might  have  caused  the  discrepancy.  First,  a 
bad  measurement  device  might  have  caused  the  lower  temperature  measure¬ 
ments,  and  this  in  turn  caused  the  lower  heat  rate.  Second,  the  non- 
isothermal  wall  effect  could  have  contributed  to  the  lower  heat  rate. 

(Ref.  4)  Since  the  leading  edge  is  stainless  steel  which  is  a  heat 
sink  type  material,  temperatures  on  the  stainless  steel  should  be  less 
than  the  temperature  on  the  HRSl  tile.  Therefore,  the  lower  temperature 
at  the  leading  edge  could  affect  the  temperature  on  the  HRSl  which  in 
this  case,  lowers  the  heat  rate  on  the  HRSl.  Nonetheless,  these  two 
possibilities  needed  further  investigation. 

The  heating  parameter  estimates  of  maneuver  one  with  the  Kalman 
filter  turned  on  were  compared  with  the  estimates  with  the  Kalman  filter 
turned  off.  Slope  and  magnitude  of  the  heating  estimates  with  the  filter 
off  were  higher  as  shown  in  Figure  0-2.  The  slope  and  magnitude  estimate 
with  the  filter  on  agree  well  with  steady  stn*-  results  as  shown  in 
Figure  D-16.  Since  the  heating  estimates  with  the  Kalman  filter  on  match 
the  heating  estimates  with  the  steady  state  cases,  the  estimates  with  the 
Kalman  filter  are  believed  to  be  better.  To  show  this,  heating  parameters 
of  maneuver  one  were  fixed  similar  to  a  simulator  run  using  the  heating 
estimates  from  Figure  D-16.  The  result  is  shovni  in  I  igure  D-17.  liie 
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result  clearly  shows  that  the  deviation  between  the  a  priori  temperature 
estimates  and  the  temperature  measurement  is  minimized,  particularly 
the  slopes  of  the  time  histories.  Since  the  slopes  of  the  temperature 
time  histories  Indicated  the  correct  heating  rate,  the  Kalman  filter 
does  optimize  the  temperature  and  the  heating  parameter  estimates. 

FRSI  Wind  Tunnel  Results  and  Discussion 

Since  there  is  only  one  thermocouple  location  available  in  the 
HRSI  tile,  data  for  the  FRSI  was  used  because  thermocouples  were  available 
in  several  locations  throughout  the  tile.  F.rst  of  all,  a  thermocouple 
at  the  same  location  as  In  HRSI  was  needed  to  confirm  the  validity  of 
the  measurement  device  in  the  HRSI.  One  thermocouple  available  was  a 
little  closer  to  the  leading  edge  of  the  test  article.  The  thermocouple 
is  denoted  as  TC  1  in  Figure  6.  A  bondline  thermocouple  was  also  useu. 
Since  the  FRSI  has  a  smaller  number  of  nodes,  an  additional  bondlinc 
thermocouple  will  not  Increase  computational  time  too  much.  The  result 
should  give  a  very  good  Indication  of  how  good  the  thermocouple  measure¬ 
ment  device  was  used  In  the  HRSI  analysis.  Second,  thermocouples  at 
two  more  locations  were  used.  (See  Figure  6  for  TC4,  TC5,  TC7  and  TCH) 

If  the  non- isothermal  wall  in  fact  caused  the  discrepancy,  the  effect 
should  be  less  as  the  distance  between  measurement  and  the  leading 
edge  Is  larger.  From  these  two  analyses,  a  better  understanding  of 
the  discrepancy  will  be  possible. 

With  the  previous  experience,  steady  state  runs  were  analyzed 
first.  The  time  history  for  the  steady  state  runs  are  shown  in  Figure 
0-18  and  Figure  D-19.  The  heating  parameters  were  estimated  at  each 
thermocouple  locations.  The  heating  parameters  estimates  at  the  fii  :t 
location  which  is  closest  to  the  leading  edge  came  to  be  about  the 
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magnitude  as  In  the  HRSI  results.  This  has  demonstrated  that  the  thermo¬ 
couple  measurement  Is  valid.  This  also  coiUi.i;...;  that  the  discrepancy 
previously  noted  was  not  caused  by  instrumentation  errors.  The  magnitude 
of  the  heat  rate  estimates  increased  as  the  location  became  further  away 
from  the  leading  edge.  Thus,  the  results  confirm  the  assumption  for 
non- isothermal  wall  effect,  that  as  the  distance  between  the  non-isotiiermal 
wall  becomes  large,  the  effect  should  be  less.  One  other  interesting 
result  was  that  the  coating  thickness  is  very  sensitive  to  the  error  iti 
the  models  even  in  the  steady  state  runs.  In  this  case,  the  error  is 
caused  by  the  non-isothermal  wall  effect.  On  the  other  hand,  the  coating 
thickness  could  be  lower  at  the  leading  edge  r.ince  the  coating  was  applied 
on  by  hand.  Therefore,  no  conclusion  will  be  drawn  from  this. 

One  transient  maneuver  run  was  analyzed  for  completeness.  The 
time  history  is  shown  in  Figure  0-20.  The  results  also  were  similar 
to  the  steady  state  results.  The  comparison  between  FRSI  Results  and 
Thin  Skin  and  Eckert  theory  is  shown  in  Figure  0-21.  Comparing  tlie  HKSl 
and  FRSI  results  at  a  similar  TC  location,  the  heating  rate  of  the  FRSI 
is  higher  than  the  heating  rate  of  the  HRSI.  The  reason  is  that  the 
specific  heat  of  the  FRSI  coating  is  higher  than  the  HRSI  coating. 

Thus,  the  FRSI  surface  temperature  will  respond  slower  than  the  HRSI 
and  will  cause  a  smaller  temperature  step  at  the  interface  between  the 
stainless  steel.  Therefore,  the  non-isothermal  wall  effect  should  be 
less. 

To  show  the  transient  error  caused  by  the  non-isothermal  wuli 
effect  as  shown  in  the  HRSI  analysis.  The  heating  and  thermal  paiu- 
meters  were  fixed  similar  to  the  HRSI  analysis.  The  results  are  ihu.-ii 
in  Figure  0-22.  The  results  have  indicated  the  transient  effect  is  i.ss 


as  the  non- isothermal  wall  effect  is  also  less.  This  confirms  that  the 
transient  error  is  caused  by  the  non-isothermal  wall  effect  also.  The 
discrepancy  is  also  investigated  by  Cappelano  (Ref.  4),  and  his  results 
are  similar.  Furthermore,  a  similar  effect  has  shown  up  in  the  flight 
thermocouple  data  analysis.  (Ref.  11)  The  analysis  was  on  the  FRSI  on 
the  Orbital  Maneuvering  System  (QMS)  pods.  The  FRSI  location  was  approxi 
mately  six  inches  downstream  of  LRSI  tiles.  A  three  second  time  skew  was 
required  to  correlate  the  flight  data,  which  is  similar  to  the  three 
second  skew  in  the  wind  tunnel  analysis. 

One  last  investigation  on  the  sample  rate  effect  has  confirmed 
the  approximation  made  using  the  scoring  assumption.  Sample  rate  has 
been  varied  from  one  sample  per  second  to  five  samples  per  second.  The 
results  of  the  parameter  estimates  are  shown  in  Table  D-II.  The  higher 
the  sample  rate  the  better  the  convergence  to  the  optimal  parameter 
estimates. 


V.  Conclusion  and  Recommendation 

Conclusion 

The  feasibility  of  the  transient  maneuver  testing  technique 
has  been  demonstrated  using  the  wind  tunnel  HRSI  and  FRSI  data.  The 
results  also  demonstrated  the  feasibility  of  the  data  reduction  method. 

The  Kalman  filter  enhanced  the  heating  parameter  estimates  even  though 
a  modelling  error  was  identified.  Without  the  Kalman  filter,  the  heating 
estimates  did  not  agree  well  with  steady  state  results.  The  biggest 
modelling  error  was  believed  to  be  caused  by  the  non-i sothermal  wall 
effect.  This  can  explain  some  of  the  questions  raised  by  the  flight 
data  analysis.  One  example  would  be  on  the  OMS.  Furthermore,  the 
transient  lag  or  the  three  second  skew  on  the  HRSI  data  analysis  also 
was  found  to  be  due  to  non- isothermal  wall  effect.  This  can  be  eliminated 
if  the  same  material  was  used  throughout  in  the  test  models.  Nonetheless, 
the  overall  results  indicated  the  transient  maneuvers  testing  teclinique 
is  feasible.  Further  development  of  the  technique  should  enhance  the 
flight  test  capability  and  also  enhance  future  wind  tunnel  testing. 

Recommendation 

Since  the  non-i sothermal  wall  effect  caused  most  of  the  modelling 
error,  a  new  test  model  with  the  same  material  properties  should  be 
used  for  future  wind  tunnel  tests.  Different  materials  with  diflerem 
material  properties  should  be  avoided  in  order  to  avoid  the  non-i •.oiliorinal 
wall  effect.  Thermocouples  at  various  locations  are  necessary  to  oi.iitu.ii 
the  assumption  of  the  non-i  sothermal  wall  effect.  For  example,  withuut. 
the  non- isothermal  wall  effect,  the  result  should  be  consistent  at  vji ious 
locations.  One  last  thing,  even  if  the  transient  test  technique  in 
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been  found  to  be  feasible  In  wind  tunnel  or  in  flight  test,  several 
steady  state  maneuver  runs  should  be  included  in  the  wind  tunnel  test. 
This  might  be  important  because  some  parameters  might  be  best  estimated 
by  steady  state  maneuvers.  For  example,  coating  thickness  was  best 
estimated  using  steady  state  runs.  This  is  due  to  the  fact  that  the 
temperature  in  the  wind  tunnel  test  does  not  reach  equilibrium  as  fast 
as  the  flight  test.  Finally,  because  the  thermal  conductivity  depends 
on  the  pressure  in  the  test  article,  the  pressure  sensors  in  the  test 
articles  should  be  sealed  to  minimize  the  uncertainty  in  the  estimates. 
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Appendix  A 


The  n-differential  equations  can  be  written 


U(t)  =  AU(t)  +  B  +  W(t) 


where  A  is  the  n  x  n  tridiagonal  matrix  with  elements. 
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The  sensitivity  equations  can  be  rearranged  in  a  more  convenient 
linear  state  space  form. 

S  =  C  S  +  D  I  '  -* 

where  C  is  a  n  x  n  tridiagonal  matrix,  and  0  is  a  n-dimensional  vectut 


With  this  form,  a  solution  of  the  sensitivity  ljii  tt.  tound  readily. 
For  computer  implementation  purposes,  finite  differences  were  u^ed 
to  approximate  the  sensitivity.  The  approximation  is 


At 


where  is  the  sensitivity  at  previous  time  t-1. 
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Table  D-2  Parameter  Estimates  for  Various  Sample  Rates 
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Figure  D-1.  Deflection  Angle  and  TeinperaLure  Tine 
History  of  Maneuver  One 
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Figure  0-14  Temperature  Time  History  of  llaniuvur  Five  wii! 

Parameters  Fixed  Similar  to  a  Simulator  Run 
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